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Three-dimensional direct numerical simulations of an incompressible open square cavity flow are 
conducted. Features of the permanent (non- linear) regime together with the linear stability analysis 
of a two-dimensional steady base flow are discussed. Spanwise boundary conditions are periodic and 
control parameters set such that the shear layer is stable against Kelvin-Helmholtz modes. Three 
branches of destabilising modes are found. The most destabilising branch is associated with steady 
modes, over a finite range of spanwise wavenumbers. The two other branches provide unsteady 
modes. Features of each branches are recovered in the permanent regime: wavelength of the most 
powerful spanwise Fourier mode, swaying phenomenon, angular frequencies, indicating that modes 
of each branches are selected and interact in the permanent flow. 



I. INTRODUCTION 

Open cavity flows belongs to the family of impinging 
flows, which are known to develop self-sustained oscilla- 
tions. When driven by a lid, cavity flows do not exhibit 
such self-sustained oscillations. Production of tones re- 
quires the existence of a shear layer, that forms between 
the outer driving flow and the inner flow, velocities in the 
inner flow being one order of magnitude smaller than in 
the outer flow, in the permanent regime. For decades, an 
extensive literature has been devoted to the study of self- 
sustained oscillations occurring in such flows in general, 
and in open cavity flows in particular. Among the most 
striking results, one may cite earlier works by Rossiter, 
who derived a phenomenological law for explaining the 
frequency selection for the the so-called Rossiter modes, 
in compressible flows |27j , Gharib and Roshko who made 
distinction between shear-layer (resonant) modes and 
wake modes [111 ]. Miksad who demonstrated the ampli- 
fication and non-linear coupling of modes excited in the 
free shear-layer [20|, HH, Rockwell and co-authors, who 
intensively studied experimentally impinging and cav- 
ity flows, and identifying several regimes of oscillations 
[151 12414261 ]. Among other features, they have explained 
the occurrence of amplitude modulations, at impinge- 
ment, in which lower frequencies are involved. Recently, 
Delprat derived an empirical formulation for explaining 
spectral features of cavity flow oscillations for various 
control parameters, in compressible regimes, that may be 
extended to the incompressible regime [5|, |6[ . Recent in- 
vestigations have considered cases of non-linear competi- 
tion of modes of oscillation, resulting in a so-called mode 



switching phenomenon, both experimentally [HI, [23j and 
numerically [7]. In recent studies, coherent structures as- 
sociated with self-sustained modes have been investigated 
in experiments by means of modern high-speed particle 
image velocimetry techniques Our present under- 
standing of open cavity flows, in the incompressible limit, 
may be roughly summarized as follows. The shear layer 
that forms between the outer and inner flows is known to 
be unstable with respect to streamwise Kelvin-Helmholtz 
modes, beyond some critical value of the Reynolds num- 
ber. Kelvin-Helmholtz modes develop, downstream in 
the shear layer, from perturbations initiated at the lead- 
ing corner. The unstable shear layer rolls up on itself to 
produce vortices that cyclically collide with the trailing 
edge. Changes in the pressure field at impingement are 
instantaneously fed back to the leading edge, due to in- 
compressibility or, equivalent ly, when the cavity length is 
small in regard to the acoustic wavelength and the mate- 
rial velocity small before sound velocity. Feedback then 
initiates the growth of a new perturbation at the lead- 
ing edge. Therefore, both leading and trailing edges are 
locked in phase, giving rise to strong self-sustained oscil- 
lations, whose frequency both depends on the incoming 
flow velocity, the cavity length, and the boundary layer 
thickness at the leading edge. 

What has been far less studied is the inner flow, 
which exhibits non-trivial three-dimensional organisation 
0, [H, 0, [22|. Centrifugal instabilities were shown to 
give rise, beyond a critical Reynolds number, to Taylor- 
Gortler-like rolls [8|, |9[. The flow curvature along the 
trailing wall bends the pathline of material particles, 
which are subject to inertial forces. When viscous terms 
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cannot overcome centrifugal effects, instabilities may de- 
velop, eventually giving rise to the formation of stable 
raw of vortical rolls, in the spanwise direction, as ob- 
served in experiments or direct numerical simulations of 
the flow [ll, 14, [9|. In experiments, as well as in three- 
dimensional direct numerical simulations of the flow, 
where spanwise walls are rigid, rolls drift toward one or 
the other wall, depending on their initial spanwise loca- 
tion [8, 9]. Drift was interpreted as initiated by a pump- 
ing effect at the walls. More recently, Alizard et at [1[ 
conducted direct numerical simulations of the cavity flow, 
with periodic spanwise boundary conditions, in a configu- 
ration where the shear layer has become unstable against 
Kelvin-Helmholtz modes. In this study, a linear stability 
analysis has been conducted, on the (unstable) steady 
base flow, with respect to spanwise Fourier modes, and 
several branches of (linearly) growing modes were found. 

The present contribution relies on the intrinsic features 
of the inner flow, in an open cavity flow at a Reynolds 
number below the threshold of Kelvin-Helmholtz insta- 
bility. Three-dimensional direct numerical simulations 
of the cavity flow are conducted with periodic spanwise 
boundary conditions, in order to get rid of the Eck- 
man pumping effect. In the case under study, though 
the shear layer remains stable against Kelvin-Helmholtz 
modes, the inner flow exhibits a three-dimensional organ- 
isation. Beside the formation of a spanwise raw of pairs of 
vortical structures, a detailed analysis of the permanent 
regime shows that the regime is not steady. Unsteadi- 
ness comes from a slightly swaying motion of the main 
structure, very similar to the one observed in A lin- 
ear stability analysis is conducted, on a two-dimensional 
steady state with respect to three-dimensional perturba- 
tions. As it will be shown, three branches of (spanwise) 
growing modes are found, whose features are recovered 
in the permanent regime with little distortion, though 
coupled through non-linear terms. 



II. CAVITY FLOW FEATURES 

Instabilities in the shear layer primarily depend on the 
dimensionless cavity length Yq = L/9o, ratio of cavity 
length, L, to boundary layer momentum thickness, Qq, 
at the leading edge. Reynolds numbers based on L and 
incoming velocity Uo (at the upstream cavity wall), is 
another relevant control parameter. For a lower part, 
cross-stream aspect ratio Tl = L/H, defined as the ratio 
of cavity length, L, to cavity height, H, is appropriate 
when considering cavity inner-flow interactions with the 
flapping shear layer. Finally, Strouhal numbers based on 
incoming velocity, Uo, and L, constitute an adequate nor- 
malisation for frequencies corresponding to shear layer 
self-sustained oscillations. 

The curvature induced by the main flow recirculation, 
inside the cavity, shown in Figure [TJi from an experimen- 
tal flow, is responsible for the development of centrifu- 
gal instabilities, over some range of the control parame- 
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FIG. 1. Instantaneous visualisation of the flow seeded with 
smoke, for Yl = 1, Ys = 6, Re# = 2 300. (Top) Cross-stream 
view in a plane (x,y), (bottom) top- view in a plane (x,z), 
located inside the cavity at y/H = —0.3. One can distinguish 
the shear-layer and the main recirculation in the inner-flow 
(top) and a (horizontal) spanwise cut of Taylor- Gortler rolls 
revealed by the seeding smoke (bottom) . 



ters. Relevant parameters, there, are the Reynolds num- 
ber based on cavity height, Ren = UoH/v, and Yl. In 
|9[, it was shown that an unambiguous parametrization 
of centrifugal instabilities requires a third dimensionless 
number; spanwise ratio S/H was chosen, although 5, the 
cavity span, mainly contribute to the selection of number 
of spanwise structures. Rather, H/0q would be a better 
choice, although not always easy to estimate in exper- 
iments. In the non- linear saturated regime, centrifugal 
instabilities give rise to vortical structures that form a 
spanwise alley of pairs of counter-rotating vortices. A 
spanwise cut of this alley of vortices, in a top view of 
an experimental flow, is shown in Figure [TJd, in a cavity 
where Yl = 1. Each pair is shaped like a torus around 
the main recirculation flow [H, Q. Spanwise walls, at 
z/S = ±1/2, generate a Bodewadt pumping, responsi- 
ble for a drift of the Taylor- Gortler rolls toward one wall 
or the other, as can be seen in the experimental space- 
time diagram of Figure [2j In this figure, oblique lines 
correspond to drifting structures, and it is clear, from 
the different slopes that can be seen, that inner spanwise 
dynamics is associated with several drift velocities, de- 
pending on the spanwise position, z. Henceforth, drift 
motions introduce additional unsteadiness in the flow, 
enriching the low-frequency range of power spectra. 
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FIG. 2. Top: spanwise topview of the flow. Bottom: space- 
time diagram of the Taylor- Gortler rolls dynamics. The space- 
time diagram is made of the vertical concatenation of the 
spatial line shown in the top picture (the line intercepts the 
mushroom- like structure, close to the cavity upstream wall). 
From an experimental cavity flow, at Yl — 1, Rch — 5 800. 



III. DIRECT NUMERICAL SIMULATIONS 

A three-dimensional direct numerical simulation is per- 
formed with the cavity aspect ratio set to Tl = 1 and 
Reynolds number to Re = 3 850. Spanwise boundary 
conditions are periodic, to get the intrinsic features of 
the inner flow. 



A. Numerical methods 



of the form: 

f/_^v 2> )*y n+1 = _vp n+1 + ^ n ' n - 1 

< V 3#e J 

y.jyn+l = q 
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where superscript n tags time t ni At is the time step and 
gn,n-i j g ^ e source term gathering all explicit quanti- 
ties, evaluated at times t n and t n -\. For each time-step, 
the numerical procedure is splitted in two parts, a pre- 
diction step and a projection step. The former consists 
in resolving the Helmholtz equation by considering the 
explicit pressure field P n in place of the implicit one. 
The integration is performed with an ADI (Alternating 
Direction Implicit) method [14[. As a result, we obtain 
an estimated velocity field U* that is not yet divergence- 
free. The incompressibility property is imposed by using 
an incremental projection method [13]. The projection 
step requires to resolve a Poisson-type equation, using 
a relaxed Gauss-Seidel method coupled to a multigrid 
method, in order to accelerate convergence, where the 
source term relies on non-zero divergence of the predicted 
velocity field: 

V 2 (j) = V • u* 

Solution <j) corresponds to the pressure time-increment, 
gradient V(j) is the correction term such that the velocity 
field is divergence free at time t n+1 . The Poisson equa- 
tion is commonly solved with Neumann-type boundary 
conditions, where the normal derivative on the domain 
limits is zero. By doing so, the boundary condition, on 
the corresponding normal velocity component, is not af- 
fected by the correction term. 



The study addresses an incompressible and isothermal 
flow whose governing equation can be described by the 
non-dimensional Navier- Stokes equations: 
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where U is the velocity field and P the pressure field. 
Numerical simulations are performed with the OLORIN 
code developed at LIMSI, which is based on an incre- 
mental prediction - projection method, see [10] for more 
details. Momentum equations are discretised with a fi- 
nite volume approach on a staggered structured grid. 
The spatial discretisation of fluxes is carried out with 
a second-order centred scheme in a conservative form 
and time derivation is approximated by a second-order 
backward differentiation formula. Viscous terms are im- 
plicitly evaluated whereas convective fluxes are explicitly 
estimated at time t n+1 by means of a linear Adams- 
Bashford extrapolation. The discretised form of the 
Navier- Stokes equations yields a Helmholtz-type problem 



B. Geometric setup 

The geometric setup consists of an open cavity capped 
with a parallelepipedic duct in which is generated the 
channel flow driving the inner cavity flow. A cartesian 
coordinate system (x, y, z), for streamwise, crosswise and 
spanwise directions, respectively, is set midspan at the 
top of the upstream cavity wall. The cavity dimensions 
are L = 5 cm, L/H = 1, S/H = 6. Upstream and 
downstream lengths of the duct are respectively L u /L = 
1 and Ld/L = 3 and its height is H v /H = 3. The total 
domain is meshed on 160 x 128 x 192 nodes, among which 
64 x 64 x 192 are devoted to the cavity. The mesh, regular 
in the spanwise direction, is particularly refined close to 
the walls and at the cavity-top in order to enhance the 
spatial resolution of boundary layers and shear-layer. 

The inlet flow is determined by Dirichlet boundary 
conditions. In order to limit the numerical domain size, 
and therefore CPU time-consumption, the upstream vein 
length has been reduced. The inlet velocity profile is 
preliminary calculated by means of a 2D simulation of 
a laminar channel flow in spatial development, represen- 
tative of the experimental upstream vein. The profile is 
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FIG. 3. Three views of Taylor- Gortler structure, as isosur- 
faces of dimensionless helicity hUq/H (edge values ±5.8). 
Only one half of the spanwise direction is represented for 
the sake of clarity. Spanwise boundary conditions are peri- 
odic. Top pictures are side-views of the structure, slightly 
seen from above and downstream (flow direction is indicated 
by the arrow). Bottom picture is a front view from the cav- 
ity downstream wall. Black corresponds to negative values of 
helicity, white to positive values. 



then extracted out of the appropriate cross-section of the 
channel-flow and extruded in the periodic spanwise direc- 
tion. The Reynolds number is set at Re = 3 850 (bulk 
velocity U q = 1.2 m/s), and the boundary layer momen- 
tum thickness is 0o/H = 26 x 10 -3 . Usual non-sliding 
conditions are applied at the walls. The numerical sim- 
ulation is carried out over a time-duration of 60 s (1440 
time-units H/U q ), after overtaking the numerical transi- 
tory state of the flow. 



C. Intrinsic features of Taylor- Gortler vortices 

Similarly to what is observed in spanwise wall-bounded 
configurations, raw of Taylor- Gortler-like vortical struc- 
tures take place inside the cavity flow and wind around 
the main recirculation vortex, as shown in Figure 03 
However, contrary to the wall-bounded configuration, no 
structure migration towards the spanwise walls is ob- 
served. 

Contrary to what could be inferred from helicity rep- 
resentation of Figure [3j where tubes of Taylor- Gortler 
vortices appear broken and the sign of helicity changes 
while moving along the tube, the material tubes are not 
broken and do form closed tori around the main recircu- 
lation (see streamlines of Figure |4]). The change of sign in 
helicity, h = lj • v, is well understood when noticing that 
helicity is made of two contributions, h = uj\\-v\\+lj±-v±, 
where uj\\ is the vorticity component along the tube axis, 





FIG. 4. Two views of streamlines stem from a spanwise line, 
at the bottom of the cavity, close to the upstream wall (the 
seeding line is materialised in pictures). Trajectories forward 
in time form inner loops (gray on top picture), trajectories 
backward in time form larger loops (black on top picture). 
Isosurfaces of Figure [3] are shown in transparency. One hun- 
dred particles seed the flow, trajectories are computed by a 
Runge-Kutta 4 time-integrator. 



#11 the velocity component parallel to aJy, and uj± and v± 
are the vorticity and velocity components, perpendicular 
to the tube axis, respectively. While cJy and m keep the 
same orientation all along the Taylor- Gortler tube, rela- 
tively to each other, uj± • v± can change of sign along the 
tube. 

Although the flow is not fully steady, it is insightful 
to consider streamlines released backward and forward 
to a seeding line at the upstream bottom quadrant of 
the cavity (materialized line in Figure 2]), and consider 
streamlines as approximate material trajectories. Tra- 
jectories gather into beams and draw funnel-like shapes, 
materializing the coherent structures highlighted by iso- 
surfaces of helicity (see Figured)). In addition, although 
particles may start aligned in the spanwise direction z, 
their trajectories coil around each other, which reveals an 
helical flow in Taylor- Got ler vortices. There is no mass 
flux between two adjacent counter-rotating vortices. The 
only spanwise transfer occurs inside both corner vortices, 
at the bottom of the cavity. Particles seeded inside the 
core of the main recirculation describe concentric circles, 
characteristic of a solid rotation (not shown in the Fig- 
ure). 

A spectral Fourier analysis on velocity components, in 
the spanwise direction, yields another insight on the spa- 
tial arrangement of the flow (see Figure [5]). The most 
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FIG. 5. Each (x, y) slice represents the power structure, in the 
(x,y) plane, of spanwise Fourier modes e tkmZ from a discrete 
Fourier decomposition of the permanent flow, where k m = 
2nm/ S is the m th mode wave- vector. Slices of fifty modes are 
shown and connected by isosurfaces with values 1.2u x /U q = 
10~ 4 (blue), 10~ 3 (yellow), 10~ 2 (red). Most energetic modes 
are found for m — 15 and m — 30 (red colors). 




X/H 



FIG. 7. Power spectral density computed from velocity time- 
recording at a probe inside the cavity. 

mainly observed on the inner blobs of helicity, close to 
the upstream cavity wall (better seen on the top right 
picture). 

It is worthwhile noticing that the mode with largest 
growth-rate found in [4], for Tl = 1 and Mach num- 
bers M > 0.3 has a wave number X/H ~ 0.5, similar 
to our X/H = 0.42. This mode, in [4(, is expected to 
be a steady mode (u = 0). In our configuration, the 
non-linearly saturated regime is found slightly unsteady, 
with ujH/U q = 0.16. This discrepancy motivates a three- 
dimensional linear stability analysis of the flow in the 
incompressible limit (M = 0). 

IV. LINEAR STABILITY ANALYSIS 

Stability properties of a two-dimensional steady base 
state, with respect to spanwise perturbations, are now 
considered. 



FIG. 6. Streamwise structure of the most energetic spanwise 
Fourier mode, m — 15 in Figure [5] Colorscale encodes power 
spectral density of u x , in units of U q . 



energetic Fourier modes, in space, are found for span- 
wise wavelengths, X/H ~ 0.42 (A/5 ~ 0.07) and 0.84 
(X/S ~ 0.14). Wavelength X/H = 0.42 is therefore as- 
sociated with the rake of vortex pairs, while its second 
harmonic, X/H = 0.84, is connected to single vortices 
inside pairs. The streamwise (x, y) structure of the most 
energetic mode is shown in Figure [6] 

In addition, a spectral Fourier analysis on temporal 
samplings of velocity, at probes set in the cavity, points 
out that the flow is not quite stationary and contains a 
slight oscillating component, as shown in Figure [7] This 
oscillation, found all over the cavity, is characterized by 
a angular frequency, uH/U q = 0.16. A flow inspection 
reveals that this oscillating mode is associated with a 
swaying motion of Taylor- Gortler-rolls, around the main 
recirculation. In Figure [3j the slightly swaying motion is 



A. Principle 

The method is described in [19]. It is aimed at charac- 
terizing the time-evolution of infinitesimal perturbations, 
u', with respect to an unstable steady base flow, L7~o, by 
means of linearised Navier- Stokes equations. The main 
instability features are depicted by the leading eigenpairs 
of the linear evolution operator, namely eigenpairs with 
the greatest real part. Leading eigenpairs are determined 
by an Arnoldi method. To shortly remind the procedure, 
consider the Navier- Stokes equations, liner arized around 
Z7o, describing the dynamics of the perturbation: 

{du' 1 
— + E7o-Vii / + u / -VE7 = -Vp' + —Au' 
ot Re 
V-u' = 0, 

(2) 

In compact form, the system rewrites: 
du' 

—-=Nuu' + Du' = Au', 
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where Njj is the linearised operator of evolution asso- 
ciated with the eulerian part of the Navier- Stokes equa- 
tions and D is the linear operator associated with viscous 
terms. Time-evolution of the perturbation is therefore 
given by: 



u 



n+l 



exip(AAt)u' n = B u' n 



(3) 



where n and n + l are the temporal subscripts referring 
to two consecutive samples with At apart and B is the 
discretised linear operator that governs the temporal evo- 
lution of u' . Operators A and B have approximatively 
the same eigenvectors while their eigenvalues are related 
by: 



Hb = exp(^Ai). 



(4) 



Therefore, leading eigenvalues \±a, of A, also provide 
dominant eigenvalues of £?, and reciprocally. Eigen- 
values and their associated eigenvectors, are deter- 
mined by considering the Krylov subspace, Uk, spanned 
by K vectors, u' k , the dimension of which is 37V, N being 
the number of grid nodes. Dimension K is small and di- 
rectly related to the number of leading eigenpairs of inter- 
est. Vectors u' k are generated by the same time-stepping 
procedure as the one used to perform the 3D DNS, with 
the usual Navier- Stokes equations (pp) replaced by their 
linearised counterpart ((3]). Time- integration is carried 
out starting from an initial condition u' that will be dis- 
cussed further. The time stepping procedure runs until 
the most decaying eigenvectors have vanished. The K 
vectors u' k are recorded at a sampling time, At, chosen 
as a multiple of the numerical time-step, At = aSt. Pa- 
rameter a is small (usually unity) and adjusted such as 
to produce results insensitive with respect to the choice 
of At. Time- integration is not conducted for too long 
times in order to avoid the solution u' k to be completely 
dominated by the most dominant eigenvector. The K 
vectors are ortho- normalised by a Graham- Schmit pro- 
cedure. The new Krylov subspace is composed of K 
orthonormalised three components vectors Vk, gathered 
into a 3N x K matrix Vk- Moreover, from the QR fac- 
torisation of Uk, it is possible to determine a relation 
yielding an approximation of BVk in the base Vk, of the 
form: 



BV K = V K H + r K+u 



(5) 



which is the Arnoldi equation. Matrix H is a K x K 
Hessenberg matrix, whose elements are the inner prod- 
ucts (vi,Bvj). A reference criterion for the selection of 
dimension K is based on the minimisation of the residual 
tk+i- The diagonalisation of H yields eigenvalues and 
eigenvectors ipk. The former approximate the dominant 
eigenvalues of B, and therefore yields leading eigenvalues 
of A from Eq. [U The later provide the corresponding 
eigenvectors <pk of B and A, given by 




FIG. 8. Two-dimensional steady base- flow, solution of equa- 
tions © for Re = 3 850 and Y L = 1. 



B. Choice of the base flow and initial perturbation 

form 



A preliminary two-dimensional numerical simulation, 
in the streamwise plane (x,y), for the same Reynolds 
number, Re = 3 850, exhibits an asymptotically stable 
steady-state flow, Uq D . As a consequence, in the three- 
dimensional configuration under study, the instability of 
the base- flow must occur with respect to spanwise modes, 
as it is observed in three-dimensional numerical simula- 
tions of the Navier- Stokes equations (pQ). Therefore, the 
two-dimensional base-flow, Uq D , is perturbated with ini- 
tial conditions of the form: 



u'(x, y, z) = uq(x, y) exp (i~^^J ex P 



(fijt) , 



that is, a mode of wavelength Xj — L/m, m G N*, 
in the spanwise direction, and (complex) growth-rate 
fij = Gj + iujj , with <jj the temporal growth-rate and ujj 
the angular frequency. Only stationary wave-like modes 
are considered since no spanwise drift is expected in the 
permanent flow. The governing equations (|3|) rewrite as: 



A(u' ) 



dt 



Nuu' 



D 



/(x,y) 



2tt\ 



0' 



where Dj^ x ^ y \ is the restriction of D to plane (x,y). For 
each specific wavelength Aj, the K leading eigenpairs 
j^j,k^j,k can be sought out. 



C. Linearly growing modes 

Linear stability analysis is carried out for a set of initial 
perturbations X/H in the range to about 3/2, corre- 
sponding to 5/4. As a first result, it is found that span- 
wise modes with vanishing wave numbers have negative 
growth-rate. Consequently, the two-dimensional steady 
base-flow is unstable with respect to spanwise, rather 
than purely streamwise, instabilities. Three branches of 
instability are found, whose temporal growth-rate a and 
angular frequency lj, as functions of A, are illustrated 



FIG. 9. Dimensionless growth-rate (top) and angular fre- 
quency (bottom) with respect to dimensionless spanwise 
wavenumber X/H. Three families are found: i) steady grow- 
ing modes (black squares), ii) oscillatory growing modes, with 
a constant angular frequency, ujH/U q ~ 0.16, and small wave- 
length, X/H < 0.62 (green triangles), Hi) oscillatory growing 
modes of larger wavelength, with an angular frequency in- 
creasing with wavelength (red circles) . Mode with the highest 
growth-rate belongs to family i) with a dimensionless wave- 
length X/H ~ 0.4. 




FIG. 10. Streamwise structure (modulus) of the spanwise 
mode with largest growth rate of the first family (oj = 0), 
referred to as mode (i). 




X/H 

FIG. 11. Streamwise structure of the spanwise eigen- function 
with largest growth-rate of the second family (ujH/U q = 
0.16), referred to as mode(ii): real part on the top, imagi- 
nary part on the bottom. The colorscale encodes u x , in units 
oiU q . 



in Figure [9j One branch (black squares) is associated 
with a stationary bifurcation (uj = 0), on the wave- 
length range X/H G [0.24,0.618]. The largest growth- 
rate on this branch is found at X/H = 0.402, close to 
the wavelength found for the raw of Taylor- Gortler vor- 
tices (Xtg/H ~ 0.42). We will refer to this mode as 
mode(i). A spatial representation of mode (i) is shown 
in Figure [TOJ 

In Figure [9] is also found, over the range of wavelength 
X/H G [0.24,0.642], another branch of solutions (green 
triangles) with positive growth-rates, and angular fre- 
quencies roughly constant over the range, at uouH/U q = 
0.16, very close to the one of Taylor- Gortler vortices. For 
this branch, the highest growth-rate is found, again, for 
X/H 0.4, and we will refer to this mode as mode (ii). A 
comparison of both real and imaginary parts of mode (ii), 
in Figure [TTJ reveals a swaying motion of the structure 
around the main recirculation. 

When the growth-rate of the stationary branch crosses 
zero, at X/H = 0.618, a new branch rises, on the 
range X/H G [0.618,1.32], characterized by a non-zero 
angular frequency (red circles). Following this oscil- 
latory branch, the angular frequency increases from 
to about uoH/U q = 0.092, with a linear variation be- 
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yond X/H = 0.78. In this family, the mode with the 
largest growth-rate, referred to as mode (iii), is found for 
Xm/H ~ 0.82 with uj Ui H/U q ~ 0.075. 

Finally, no spanwise mode is found with positive 
growth-rate when X/H > 1.32. 

D. From instability to permanent flow 

The inner flow most energetic structure, in the perma- 
nent flow, has both a spanwise wavelength, Xtg/H, and 
a streamwise structure, in close accordance with those 
of mode(i), the most linearly unstable mode in our lin- 
ear stability analysis (compare Figures l6l and [T0|) . How- 
ever, mode (i) is steady, whereas the permanent regime 
exhibits a slightly unsteady swaying motion. Note that 
transient dynamics, although not fully excluded, are very 
unlikely, since simulations in the permanent flow cover a 
time-range of 1440 time-units H/U q . In fact, when con- 
sidering the power spectral density of Figure [7J in the 
permanent regime, a main peak is found at UJTGH/U q ~ 
0.16, which precisely is the angular frequency of modes of 
the second family. Moreover, the streamwise structure of 
mode (ii), in that family, exhibits a swaying-like motion, 
as shown in Figure [TTJ when comparing both real and 
imaginary parts. Such features suggest that this mode, 
as mode (i), is also selected by the flow. Note that linear 
modes are recovered in the permanent regime without 
significant distortion, despite of non-linear effects. 

It is worth to note that the angular frequency of 
the third family mode, with the highest growth-rate, is 
uJiuH/Uq ~ 0.07. Would such a mode be present in the 
flow, it could couple to mode (ii) and generate compo- 
nents at (ujh ± 0Jm)H '/Uq, i.e. 0.09 and 0.23. In the 
power spectral density of Figure [7| a lobe is actually 
found around ooH/U q ~ 0.25, suggesting a non-linear 
(quadratic) coupling of mode (iii) with mode (ii). In addi- 
tion, the spanwise Fourier mode associated with m = 30, 
in Figure [5] for the permanent regime, is not found in 
the linear stability analysis, indict aing that this mode is 
rather generated by quadratic non-linear self-interaction 
of mode m = 15. 



ible limit and at Reynolds number Re = 3 850, below the 
threshold of shear layer instability. They reproduce most 
of the inner flow characteristic features observed in ex- 
periments. As a result of centrifugal instabilities, vorti- 
cal structures develop and form spanwise alley of pairs of 
counter-rotating vortices [§|, [9j . When spanwise bound- 
ary conditions are periodic, the alley of vortices is fixed. 
However, in the non-linear ly saturated regime, a closer 
inspection of the flow reveals a slightly swaying motion 
of the structure of vortices around the main recirculation 
flow. 

The fastest (linearly) growing spanwise mode, in a lin- 
ear stability analysis of the steady base state, has a wave- 
length X/H = 0.4, fairly close to the wavelength associ- 
ated with the spanwise alley of vortices. However, from 
stability analysis, this mode is expected stationary. Yet, 
the swaying motion is caught by a mode of the second 
family, whose growth-rate is maximal for wavelengths of 
the order of Xtg- Surprisingly, streamwise structures, 
wavelength and angular frequencies of linear modes are 
not much distorted in the non-linear regime. Non-linear 
effects are revealed by coupling modes of two families, 
namely the modes with largest growth-rates from fami- 
lies (ii) and (iii). 

Note that beyond the drift motion initiated by the 
pumping effect at the spanwise walls, a secondary bi- 
furcation towards an unsteady oscillating phenomenon 
of the alley of vortices has been reported in experiments, 
see for instance [9j. Although the streakline features of 
this secondary instability could possibly be understood 
as the result of a swaying motion, it is not possible, yet, 
to ascertain whether both phenomenons, in simulations 
and experiments, be related or not. 
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